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PTEN is a tumour suppressor frequently mutated in many types of cancers. Here we show that 
targeted disruption of PTEN leads to neoplastic transformation of human neural stem cells 
(NSCs), but not mesenchymal stem cells. PTEN-deficient NSCs display neoplasm-associated 
metabolic and gene expression profiles and generate intracranial tumours in immunodeficient 
mice. PTEN is localized to the nucleus in NSCs, binds to the PAX7 promoter through association 
with cAMP responsive element binding protein 1 (CREB)/CREB binding protein (CBP) and 
inhibits PAX7 transcription. PTEN deficiency leads to the upregulation of PAX7, which in turn 
promotes oncogenic transformation of NSCs and instates 'aggressiveness' in human glio- 
blastoma stem cells. In a large clinical database, we find increased PAX7 levels in PTEN- 
deficient glioblastoma. Furthermore, we identify that mitomycin C selectively triggers apoptosis 
in NSCs with PTEN deficiency. Together, we uncover a potential mechanism of how PTEN 


safeguards NSCs, and establish a cellular platform to identify factors involved in NSC trans- 
formation, potentially permitting personalized treatment of glioblastoma. 
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suppressor whose loss-of-function mutations are often 

encountered in human cancers. PTEN mutations are 
observed in 60% of glioblastoma multiforme (GBM) and are 
among the most frequent genetic alterations linked to GBM!. 
GBMs bearing PTEN loss-of-function mutations are usually 
associated with increased invasive behaviours and drug 
resistance~*. Glioblastoma stem cells (GSCs), the tumorigenic 
component of GBM, represent a rare cell population that are 
resistant to conventional radio- or chemo-therapy, and are 
presumably involved in cancer relapse?$. Evidence from mouse 
tumour models reveals neural precursor/stem cells as the cell-of- 
origins for GBM or GSCs’-’, and GBM is postulated to be 
derived from transformed neural stem cells (NSCs) that undergo 
carcinogenic hits!®. The high mutation rate of PTEN in GBM 
suggests its potential as one of the initiating oncogenic events or a 
key factor in promoting cancer aggressiveness, similarly as seen in 
endometrial cancer!!. The correlation between PTEN deficiency 
and poor prognosis suggests a more complex role of PTEN loss in 
GBM progression. These observations raise an interesting 
question, that is, how PTEN loss leads to GBM initiation or 
promotes its progression? 

Mouse models have been successfully used to investigate the 
roles of genetic mutations in triggering oncogenic NSC 
transformation and/or mediating GBM pathogenesis!^P?, The 
known differences between mouse and human cancer biology, 
including differential telomere length, distinct utilization of 
pl6lNK4 RB versus p53 signalling and different sensitivity to 
anti-tumour drugs, however, have limited the degree to which 
insights derived from mouse models can be directly translated 
to human applications!^-16, The advances in human stem cells 
and targeted gene editing technology have opened a new avenue 
for disease modelling and drug discovery". Although many 
genetic disease models that are linked to development and ageing 
have been developed using human embryonic stem cells (ESCs) 
or induced pluripotent stem cells (iPSCs)!/ 74, very few human 
cancer models employing targeted genetic mutations in adult 
stem cells have been established for gaining mechanistic insights 
or testing drug efficacies^?26, 

Considering the potential of NSCs being the cell-of-origin 
for human GBM, and PTEN deletion has frequently been 
reported in GBM, we hypothesize that PTEN functions as 
a gatekeeper to protect human NSCs from neoplastic trans- 
formation. Accordingly, we generated PTEN-deficient human 
NSCs by targeted gene editing. PTEN deficiency resulted in 
a reprogramming of NSCs towards a GSC-like phenotype in 
a highly lineage-specific mechanism primarily through 
transcriptional activation of PAX7, a pathway that is also linked 
to the aggressive characteristics of GSCs in GBM patients. 


p hosphatase and tensin homolog (PTEN) is a potent tumour 


Results 

PTEN deficiency leads to neoplastic features in NSCs. To 
generate isogenic PTEN-null human NSC lines, we first utilized 
transcription activator-like effector nuclease (TALEN)-mediated 
homologous recombination (HR) to delete the exon 1 of PTEN 
gene in human ESCs (Fig. 1a). Successful gene targeting at PTEN 
locus was verified by genomic PCR (Fig. 1b) Immuno- 
fluorescence staining revealed a punctate staining pattern of 
PTEN in the nucleus of wild-type (WT) ESCs, which was absent 
in PTEN homozygous knockout (PTEN -!-) 
ESCs (Fig. lc). Western blot further confirmed the loss of 
PTEN protein in PTEN-/- ESCs (Fig. 1d). The deficiency of 
PTEN, however, did not alter ESC-associated features such as 
morphology, DNA hypo-methylation at OCT4 promoter and the 
expression of pluripotency markers OCT4, SOX2, NANOG and 


TRA-1-81 (Supplementary Fig. la-c). PTEN~/~ ESCs also 
maintained the capability of differentiating into endoderm, 
mesoderm and ectoderm lineages in vivo, as did WT ESCs 
(Supplementary Fig. 1d). 

To study the role of PTEN in NSCs, we derived NSCs from both 
WT and PTEN-/- ESCs. We confirmed the loss of PTEN in 
PTEN-!- NSCs (Fig. lef and Supplementary Fig. 2a). Both the 
WT and PTEN’ NSCs exhibited typical features of neural 
progenitors, including the expression of neural progenitor markers 
PAX6, SOX2 and Nestin, the absence of pluripotency-related 
antigen TRA-1-81, the hyper-methylation of OCT4 promoter, as 
well as the enrichment of H3K4me3 levels at PAX6 and Nestin loci 
(Fig. 1g and Supplementary Fig. 2b-f). More importantly both WT 
and PTEN~/~ NSCs were able to efficiently differentiate into 
neurons and astrocytes in vitro (Fig. 1g and Supplementary 
Fig. 2g), confirming their NSC identity. 

Next, we investigated whether PTEN-deficient NSCs could 
have acquired neoplastic potentials. PTENT’ NSCs exhibited 
accelerated growth under suspended neurosphere culturing 
condition (Supplementary Fig. 3a), accompanied by a mild 
increase in the population of cells retained in the S and G2/M 
phases and a concomitant decrease in the population of cells in 
the G0/G1 phases (Supplementary Fig. 3b). Consistent with an 
acquired growth advantage, PTEN^ /^ NSCs exhibited enlarged 
nuclei size, increased number of nucleoli and elevated rRNA 
transcripts (Supplementary Fig. 3c-e). It has been reported that 
PTEN-loss is associated with genomic instability in mouse 
embryonic fibroblasts (MEFs)7. We performed karyotyping 
and genome-wide copy number variation (CNV) analyses by 
deep sequencing and did not observe any discernible 
chromosomal aberrations in PTEN~/~ NSCs compared with 
the WT NSCs (Supplementary Fig. 3f, g). In agreement with a 
neoplastic potential, PTEN~/~ NSCs showed resistance to 
stresses-induced cell death (Supplementary Fig. 3hj)?9??, as 
well as significant elevation in cell migration and clonal expansion 
abilities (Fig. 2a,b and Supplementary Fig. 4a,b). Re-introduction 
of PTEN into the PTEN-/- NSCs partially rescued these 
neoplastic phenotypes (Fig. 2c,d and Supplementary Fig. 4c,d). In 
addition, knockdown of PTEN using small hairpin (sh) RNA in 
either ESC- or iPSC-derived NSCs recapitulated the aggressive 
phenotypes observed in PTEN-/- NSCs (Supplementary 
Figs 4ef and 5a-e). Altogether, these im vitro experiments 
demonstrated that PTEN deficiency endowed NSCs with 
neoplastic potential. 

Next we explored whether PTEN deficiency is linked to tumour 
initiation in NSCs using an in vivo mouse model. We implanted 
WT or PTEN-deficient NSCs expressing luciferase into the brains 
of immunocompromised NOD/SCID mice. PTEN-deficient NSCs 
were able to grow efficiently and form intracranial tumours 
in vivo as evidenced by both positive luminescence and magnetic 
resonance imaging (MRI) signals (Fig. 2e and Supplementary 
Fig. 6a-c). Similar to primary human NSCs*°, WT ESC-derived 
NSCs did not expand in the brain (Fig. 2e). Haematoxylin and 
eosin (H&E) staining of brain slices showed a substantial 
expansion of the grafted PTEN~/~ NSCs, with characteristics 
of neoplastic cells intruding adjacent normal brain tissue and 
histological resemblance of intracranial tumour derived from the 
grafted human GSCs (Fig. 2f and Supplementary Fig. 6e)?0?, 
Occasionally we observed invasion of neoplastic cells into other 
brain regions, intra-tumour vasculatures and non-classic necrosis 
in the mouse brain implanted with PTENT’ NSCs for a longer 
course (~70 days; Supplementary Fig. 6e). PTEN~/~ NSC- 
derived neoplasm also expressed typical mixed markers of 
GBM including Nestin, SRY box 2 (SOX2), glial fibrillary 
acidic protein (GFAP), neuron-specific beta -III tubulin (Tujl) 
and microtubule-associated protein 2 (MAP2) (Fig. 2g 


2 NATURE COMMUNICATIONS | 6:10068 | DOI: 10.1038/ncomms10068 | www.nature.com/naturecommunications 


NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10068 


ARTICLE 


Exons 1 2 345 6789 
PTEN locus + —— NNF 
a P5 
— eaa 


P1 (1,110 bp) A 1,087 bp 
PTEN locus “ATG = 
315b] P6 
>< Pe OO 
HR-donor vector 5'Homology arm, S S co Homology arm 
— 
1,185 bp P2 1,087 bp 


AExon1 (795 bp) JHR 


After targeting 
AExon1 (795 bp)| Neo removal 


After Neo removal ——J-————— — 2 ——— —— ————— 


d ESC 


PTEN** PTEN 
kDa 


50 — 
50 — 
50 — 


SOX2/PAX6/DNA 


PTEN*^ 


PTEN^ 


Nestin/Ki67/DNA 


b c 
BEN, ee x ESC PTEN/DNA 

2 2 
W W 

E E i 

" i 

z 

. [m] 

PTEN/DNA f 


NSC 


PTEN** PTEN^ 


kDa P 
PTEN 


p-PTEN 


50 
50 
50 


Neuron 


Tuj1/MAP2/DNA 


Figure 1 | Generation and characterization of PTEN-deficient NSCs. (a) Schematic representation of TALEN-based PTEN targeting strategy. Primers used 
for b are shown as arrows (P1-P6). The donor vector includes a neomycin-resistance cassette (Neo) allowing for positive selection. (b) PCR analysis of WT 
and PTEN~/~ ESCs using primer pairs indicated (P1+ P2: 1,744 bp; P3 + P4: 1,872 bp; P5 + P6: 1,035 bp). (c) Immunofluorescence analysis performed on 
WT and PTEN -/ — ESCs with an anti-PTEN antibody. PTEN was absent in the PTEN -/ ESCs. Nuclei were stained with Hoechst 33342. Scale bars, 
12.5 uum. (d) Immunoblotting verified the absence of PTEN protein in PTEN^/ — ESCs with anti-PTEN and anti-phospho-PTEN (Ser380) antibodies. 
B-Tubulin was used as loading control. (e) Immunostaining of PTEN in WT and PTEN~/~ NSCs with an anti-PTEN antibody. Scale bars, 10 um. 

(f) Immunoblotting analysis of PTEN and phospho-PTEN in WT and PTEN~/~ NSCs. (g) Immunostaining of neural progenitor- (left) and neuron- (right) 
specific markers in WT and PTEN~/~ NSCs (left) and their neuronal derivatives (right). Scale bars, 25 um (NSC) and 50 um (neuron). 


In agreement with several studies using GSCs?96, the 
neoplastic phenotype in PTEN~/~ NSCs was mediated by 
their CD133-positive, rather than CD133-negative population 
(Supplementary Fig. 6f). The in vivo aggressiveness of PTEN -/ - 
NSCs was able to be suppressed by expressing exogenous 
PTEN (Supplementary Fig. 6d), suggesting a direct 
involvement of PTEN in suppressing the neoplastic phenotype 
in NSCs. 

A full-fledged cancer phenotype normally needs a combination 
of multiple oncogenic mutations. To test whether the neoplastic 
phenotype in PTEN-deficient NSCs can be aggravated by 
inactivation of another tumour suppressor, we knocked down 
one of the most well-known tumour suppressor p53 in 
PTEN-deficient NSCs?”. As expected, we observed dramatic 
neoplastic transformation in NSCs when p53 and PTEN were 
simultaneously inactivated (Supplementary Fig. 7a-e). Next, we 
examined whether the intracranial niche of the brain was 
required for tumorigenesis mediated by PTEN deletion in 


NSCs. When PTEN~/~ NSCs were implanted into NOD/SCID 
mice via the subcutaneous route, we did not detect any neoplasm- 
like growth, suggesting that PTEN-/- NSCs required a 
supportive brain environment for tumour initiation. 

We also investigated whether PTEN deletion had impact on 
neoplastic transformation in other human adult stem cells. 
To this end, we differentiated PTENT’ ESCs to mesenchymal 
stem cells (MSCs)2%23, The PTEN-/- MSCs expressed 
specific markers, such as CD73, CD90, CD105, and were 
absent of MSC-irrelevant markers, including CD34, CD43 and 
CD45 (Supplementary Fig. 8a,b). The PTENT’ MSCs exhibited 
typical features of premature senescence, including the 
positive SA-B-Gal staining, expression of senescence-related 
gene transcripts and the compromised migration and clonal 
expansion (Supplementary Fig. 8c-g). Likewise, knockdown of 
PTEN using shRNA led to senescence in WT human 
MSCs (Supplementary Fig. 8h,i). Similar accelerated ageing 
phenotypes were also observed in the fibroblast-like cells 


NATURE COMMUNICATIONS | 6:10068 | DOI: 10.1038/ncomms10068 | www.nature.com/naturecommunications 3 


ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10068 


o 
[oz 


+ ES 8 
+ i $ kkk 
2 9 z 2 — 
i $a E B. 6 
a oS a $9 
ge ^ EE 4 
= Oo'o 
: 39 AE 
+ [ra + 8 
T l © 
Z 2 E 
z E m. 
al NSC a NSC: x oy 
OS 
ws 
c Lenti-Luc Lenti-PTEN E Lenti-Luc d Lenti-Luc Lenti-PTEN EE Lenti-Luc 
84 IB] Lenti-PTEN p m EE] Lenti-PTEN 
zh ES L3 dock 
t doe * m= 
z P E 
u o 6 ror ģ z 3 
a 2 5 " 7 ŞS 
r oO 
oS Sol 
EE os 
— 4 wt 
DI GM 9 2 T > 
" “gh $8 " yw. : E = ; 
zi * E oj ws z $T NS 
E r— u E: 
à E 
- 0 
NSC: PTEN** PTEN-- NSC: PTEN** PTEN” 
e PTEN** PTEN^ Ctrl PTEN"* PTEN^- 
f PTEN** PTEN” 
g SOX2 GFAP Tuj1 MAP2 Nestin 
Nestin Human nuclei Human nuclei Human nuclei Human nuclei GFAP 


DNA DNA DNA DNA DNA DNA 


Figure 2 | PTEN-deficient NSCs demonstrated neoplastic features in vitro and in vivo. (a) Migration abilities of WT and PTEN -/- NSCs were evaluated 
by Transwell assays. Relative cell migration efficiency was determined. Data are shown as mean t s.e.m. n= 3. ***P<0.001 (t-test). Scale bars, 1mm. 
(b) Clonal expansion analysis in WT and PTEN-/ — NSCs. Crystal violet staining-positive cells were calculated and presented as fold induction using 
Image J. Data are shown as mean € s.e.m. n= 4. ***P « 0.001 (t-test). (c) Cell migration analyses of NSCs transduced with a PTEN expression vector 
(Lenti-PTEN) or a control vector (Lenti-Luc). Data are shown as mean t s.e.m. n= 6. ***P < 0.001 (t-test); NS, not significant. Scale bars, 1mm. (d) Clonal 
expansion analysis of NSCs transduced with Lenti-PTEN or Lenti-Luc. Data are shown as mean * s.e.m. n= 3. ***P « 0.001 (t-test); NS, not significant. 
(e) Representative photon flux images from the brain of NOD/SCID mice implanted with WT (left) and PTEN~/~ (right) NSCs expressing luciferase (also 
see Supplementary Fig.6a for details). Images were taken 20 min after i-luciferin was injected intraperitoneally (i.p.). n — 3. (f) H&E staining and 
photomicrograph of a section of entire brain showing dramatic neoplastic expansion of PTEN~/~ NSCs relative to WT control (2 x 106 cells per injection), 
35 days after cells were injected into the contralateral corpus striatum. H represents residual human NSC graft without expansion, N represents neoplasm 
and P represents mouse brain parenchyma. Scale bars, 500 um. (g) Immunostaining of brain sections described in f with antibodies against the indicated 
antigens. Human nuclei, human nucleus antigen. DNA was stained with Hoechst 33342. Scale bars, 25 um. 
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derived from in vivo differentiated PTEN~'~ ESCs p53 in MSCs led to strengthened proliferation and migration 
(Supplementary Fig. 8j). The pro-senescence effect observed in instead (Supplementary Figs 7e and 8k). Herein, these 
MSCs appears to be specific to PTEN deficiency, as knockdown of observations indicate that PTEN deficiency preferentially 
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Figure 3 | Transcriptome and epigenome analyses in PTEN-deficient stem cells. (a) Glycolysis flux analysis measured by extracellular acidification rates 
(ECARs) in PTEN */ * and PTEN~/~ NSCs. n=5. (b) GO biological process analysis of genes upregulated (q-value < 0.05, FC[PTEN ^7 - /PTEN +/ 
+]> 2) after PTEN knockout in NSCs and genes downregulated (q-value < 0.05, FC[PTEN —/ - /PTEN */ * ] « 0.5) after PTEN knockout in MSCs. (c) Heat 
map displaying the upregulated genes (q-value < 0.05, FC[PTEN ^ / - /PTEN*+/+]>2) in PTEN-deficient NSCs (left) and the downregulated genes 
(q-value «0.05, FC[PTEN -/ - /PTEN */ * ] « 0.5) in PTEN-deficient MSCs (right). The representative enriched GO terms of biological processes were 
shown at the left side of each panel. (d) Boxplots showing H3K4me3, H3K27me3, 5mC and 5hmC levels in the region of 1kb upstream and 1kb 
downstream of TSS of differentially expressed genes between PTEN-/- and WT NSCs including the upregulated genes (left panels, q-value < 0.05, 
FC[PTEN -/ - /PTEN */ * ] 2 2; n— 326) and downregulated genes (right panels, q-value «0.05, FC[PTEN -/ - /PTEN */ * 1 « O.5; n — 281). *P value 

« 0.05 (Wilcoxon signed-rank test); **P value «0.01 (Wilcoxon signed-rank test); NS, not significant. 
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induces neoplastic transformation in NSCs by promoting 
proliferation and invasiveness; however, in the case of 
mesodermal cell types such as MSCs and fibroblasts, PTEN 
deletion primarily triggers cellular senescence. 


PTEN’ NSCs display oncogenic gene expression changes. 
As PTEN is a suppressor of PI3K-AKT pathway primarily via its 
activity as a lipid phosphatase????, we next examined the 
downstream signalling pathways that are regulated by PTEN. 
PTENT’ NSCs exhibited an activated AKT pathway as 
evidenced by the phosphorylation of AKT and its downstream 
targets such as mammalian target of rapamycin (mTOR) and 
glycogen synthase kinase 3 beta (GSK3B) when compared with 
the WT NSCs (Supplementary Fig. 9a,b). In contrast, AKT 
pathway activities were not affected by PTEN status in MSCs 
(Supplementary Fig. 9a), suggesting a different regulatory 
mechanism for AKT activation in MSCs. 

Given that the AKT and mTOR activations are important for 
metabolic regulation of cancer cells, we further investigated the 
metabolic changes conferred by PTEN deficiency in NSCs. Mass 
spectrometry-based screening identified 16 metabolites were 
significantly upregulated in PTEN-/- NSCs compared with 
WT NSCs (Supplementary Fig. 9c,d and Supplementary Data 1). 
Among the 16 metabolites some belong to glycolytic pathway 
including glucose, glucose-6-phosphate and fructose-6-phos- 
phate. Seahorse analysis indicated an increased glycolysis activity 
in PTEN-/- NSCs (Fig. 3a and Supplementary Fig. 9e). 
Consistently, PTENT’ NSCs produced threefold more lactate 
acid content than that produced by WT NSCs (Supplementary 
Fig. 9f), and the ratio between ADP or AMP and ATP were 
increased in PTE. NSCs relative to that of WT NSCs 
(Supplementary Fig. 9c and Supplementary Data 1). These data 
strongly support that PTENT’ NSCs have dysregulated glucose 
metabolism and favour aerobic glycolysis, a feature characteristic 
of GBM and most cancer cells*®*!, The PTEN-deficiency- 
induced metabolic change was specific to NSCs, because 
these metabolites did not increase in PTEN-/- ESCs when 
compared with WT ESCs (Supplementary Figs 9d,f and 
Supplementary Data 1). 

Next, we investigated the impact of PTEN loss on global 
transcriptomic landscape of NSCs. We performed RNA-seq 
analyses for both WT and PTEN~/~ NSCs as well as MSCs. 
Using a threshold of twofold change with q-values less than 0.05, 
we identified 607 genes differentially expressed between WT and 
PTENT’ NSCs, among which 326 were upregulated 
(Supplementary Fig. 10a-e and Supplementary Data 2). Gene 
ontology (GO) terms enriched in upregulated genes are linked to 
GSC-related functions (Fig. 3b and Supplementary Data 2). 


Notably, most of these genes belong to signalling pathways 
important for cell motion, migration and proliferation (Fig. 3b,c, 
Supplementary Fig. 10fg and Supplementary Data 2). Consistent 
with the metabolic changes, a number of enzymes involved in 
glycolysis were upregulated in PTEN~/~ NSCs (Fig. 3c and 
Supplementary Fig. 10f,g). These gene expression changes were 
partially rescued by re-introduction of PTEN in PTEN '~ NSCs 
(Supplementary Fig. 10h). In contrast, PTEN deletion in MSCs 
resulted in differential expression of 1,854 genes (Supplementary 
Fig. 10a-e). Consistent with the non-neoplastic phenotype 
in PTENT’ MSCs, PTEN deficiency led to a coordinated 
downregulation of genes required for cell cycle process (Fig. 3b,c 
and Supplementary Data 2), which is in agreement with a 
molecular signature of accelerated senescent MSCs”. Of note is 
that only few gene expression changes caused by PTEN deficiency 
overlap between NSCs and MSCs (Supplementary Fig. 10d and 
Supplementary Data 2), serving a possible explanation for the 
dramatically different cellular phenotypes observed between 
PTEN-deficient NSCs and MSCs. 

To explore whether global transcriptome change in PTEN~/~ 
NSCs are linked to epigenetic alterations, we performed 
chromatin immunoprecipitation sequencing (ChIP-seq) to map 
genome-wide distribution of H3K4me3, H3K27me3, as well as 
reduced representation bisulfite sequencing (RRBS) and 
Tet-assisted RRBS (TA-RRBS) for 5-methylcytosine (5mC) and 
5-hydroxymethylcytosine (bhmC) in WT or PTEN~/~ NSCs 
(Fig. 3d and Supplementary Fig. 11a,b). Notably, the upregulated 
genes showed significant increase of H3K4me3 and decrease of 
H3K27me3, whereas the downregulated genes were associated 
with diminished H3K4me3 and elevated H3K27me3, at their 
promoter regions (Fig. 3d). We did not observe any global 
associations between 5mC (or 5hmC) and gene expression 
changes (Fig. 3d). These observations suggested that changes in 
histone modifications including H3K4me3 and H3K27me3 could 
account for PTEN-loss-induced transcriptional dysregulation in 
PTEN~'~ NSCs. 


PAX7 mediates oncogenic features in PTENT’ NSCs. Among 
the upregulated genes, transcription factor PAX7 was at the top of 
the list and expressed 155-fold higher in PTEN~/~ NSCs than in 
WT NSCs (Supplementary Data 2). Also it is known that elevated 
levels of PAX7 in embryonal rhabdomyosarcoma cells are asso- 
ciated with increased migration and invasiveness***3. Thus, we 
focused on a potential role of PAX7 in regulating neoplastic 
transformation of NSCs. 

We confirmed the elevated PAX7 expression in PTEN~/~ 
NSCs relative to WT NSCs, using reverse transcription- 
quantitative PCR (RT-qPCR) and western blotting (Fig. 4a and 


Figure 4 | Transcriptional activation of PAX7 underlies oncogenic features in PTEN-deficient NSCs. (a) RT-qPCR analyses of PAX7 expression in WT and 
PTEN~/~ NSCs. n=3. ***P<0.001 (t-test). (b) RT-gPCR showing reduced transcripts of PAX7 in PTEN~’~ NSCs transduced with lentivirus encoding 
PTEN. n=3. ***P « 0.001 (t-test). (c) Cell migration analysis in WT and PTEN~/’~ NSCs transduced with control shRNA or PAX7 shRNA. n=5. *P« 0.05 
(t-test); and ***P< 0.001 (t-test). (d) MRI analysis of intracranially implanted WT and PTEN~/~ NSCs pre-transduced with control or PAX7 shRNA. Relative 
volumes of tumours are presented. n= 4. **P« 0.01 (t-test), and ***P< 0.001 (t-test). (e) Cell migration analysis in PTEN~/’~ NSCs transduced with the 
lentiviral vector encoding dominant negative mutants of PAX7 (DN1: a.a.1-263; DN2: a.a.1-318) or a luciferase control (Luc). n= 3. ***P<0.001 (t-test). 
(f) Cell migration analysis in WT NSCs transduced with the lentiviral vector encoding PAX7. n = 3. ***P < 0.001 (t-test). (g) ChIP-PCR showing the association 
of endogenous PTEN with PAX7 promoter in WT NSCs. The PAX7 promoter was amplified by PCR from either genomic DNA as input (lanes 1 and 4) or anti- 
PTEN immunoprecipitated DNA (lanes 3 and 6). (h) Luciferase reporter assay showed that the PTEN-docking DNA element described in g has a higher cis- 
activation ability when PTEN was depleted in NSCs. n —3. **P < 0.01 (t-test). (i) Immunoblotting analysis of CREB and phospho-CREB (Ser133) expression in 
WT and PTEN~/~ NSCs. B-Actin was used as loading control. (j) Upper panel: Identification of CRE-containing genes in the human genome by the FIMO 
software tool. Lower panel: Bioinformatic analysis predicated that 167 (51.296) of 326 genes upregulated in PTEN ^ / — NSCs are potential CREB-targets. (k) 
ChIP-gPCR analysis for PTEN, p-CREB and CBP enrichment within the PTEN-docking DNA region at the PAX7 promoter in WT and PTEN~/~ NSCs. P1-P8 
indicate the different sub-regions. ***P « 0.001 (t-test); NS, not significant. (I) Cell migration analysis (top panels) of PTEN~/~ NSCs transduced with the 
lentiviral vector encoding PTEN, PTEN(G129E), PTEN(Y138L) or a luciferase control (Luc), and the expression of PTEN was determined by 
immunofluorescence (lower panels). Scale bars, 1mm (migration assay) and 50 um (IF). n — 3. ***P« 0.001 (t-test). 
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Supplementary Fig. 12a). In contrast, the expression of PAX7 in 
PTENT’ MSCs did not increase when compared with WT 


MSCs (Supplementary Fig. 12b). The upregulation of PAX7 in 


PTENT’ NSCs was directly caused by the absence of PTEN as 
reintroduction of PTEN in these cells repressed PAX7 expression 
(Fig. 4b and Supplementary Fig. 12c). 
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To investigate whether PAX7 mediates the neoplastic effect of 
PTEN deficiency during the NSC transformation, PTEN~!~ 
NSCs were transduced with lentiviral particles encoding two 
independent PAX7 shRNAs, which resulted in downregulation of 
PAX7 in the PTEN-/- NSCs (Supplementary Fig. 12d). 
Silencing PAX7 effectively inhibited the clonal expansion and 
migratory capabilities of PTEN~'~ NSCs in vitro, as well as 
the intracranial neoplasm formation in vivo (Fig. 4c,d and 
Supplementary Fig. 12d-h). Ectopic expression of PAX7 
dominant-negative mutants, which do not contain the 
transactivation domain, also inhibited the neoplastic features of 
PTEN~'~ NSCs (Fig. 4e), suggesting a critical role of the 
transcriptional activity of PAX7 in promoting neoplastic 
phenotypes in NSCs. To examine if elevated level of PAX7 alone 
is able to cause the neoplastic changes in WT NSCs, we 
overexpressed PAX7 in WT NSCs, to a similar level as in 
PTEN~'~ NSCs, and found that ectopic PAX7 expression 
recapitulated the major phenotypes caused by PTEN depletion, 
including an increase of in vitro cellular migration and in vivo 
neoplastic expansion (Fig. 4f and Supplementary Fig. 13a-c). 
In addition, overexpression of PAX7 in WT NSCs resulted in a 
gene expression signature resembling that of PTENT’ NSCs 
and distinct from WT NSCs (Supplementary Fig. 13d,e and 
Supplementary Data 3). We found that majority of genes (19/23) 
whose transcripts were co-upregulated by PAX7 overexpression 
and PTEN-knockout are associated with tumour-promoting or 
increased cellular mobility (Supplementary Fig. 13f and 
Supplementary Data 3). Therefore, PAX7 appears to be one of 
the key factors mediating neoplastic potential of PTEN-deficient 
NSCs. 

Given that PTEN was predominantly localized to the nucleus 
in NSCs (Fig. le), and emerging evidence suggests that 
PTEN associates with chromatin DNA and is involved gene 
transcription regulation?”4, we thus reason that PTEN may 
repress PAX7 transcription through direct binding to its 
promoter. A ChIP-based promoter mapping revealed a 981-bp 
DNA region upstream of the transcription starting site (TSS) of 
PAX7 as being a docking site for PTEN in WT NSCs, but not 
in WT MSCs (Fig. 4g and Supplementary Fig. 14a). When 
cloned in a luciferase reporter, this region showed higher 
cis-activation ability when PTEN was depleted in NSCs 
(Fig. 4h). Bioinformatics analysis predicted several binding sites 
in this region for transcription factor CREB and its co-activator 
CBP/p300  (http://www.cbrc.jp/research/db/TFSEARCH.html). 
Given that CREB is a direct phosphatase substrate of 
nuclear PTEN?, we hypothesized that interplay between PTEN 
and CREB account for the PTEN-mediated transcriptional 
repression of PAX7 as well as other genes in NSCs. Supporting 
this hypothesis, PTEN-/- NSCs demonstrated elevated 
levels in the phosphorylated form of CREB (p-CREB; Fig. 4i). 
Bioinformatics predicted that 167 (51.2%) of 326 genes that 
were upregulated in PTEN~/~ NSCs are potentially CREB 
targets (P<10E~ 5; Fig. 4j). Among these genes, the promoters 
of 108 genes showed clear enrichments of the active histone mark 
H3K4me3, suggesting a permissive/active gene transcription 
state (Fig. 4j, Supplementary Fig. 14b and Supplementary 
Data 4). As for PAX7 gene, ChIP-quantitative-PCR showed 
that in WT NSCs PTEN and CREB/CBP co-occupied within 
the 981-bp DNA region at the PAX7 promoter (Fig. 4k). 
Co-immunoprecipitation (IP) and ChIP analyses further 
indicated that PTEN formed a protein complex with CREB 
and CBP (Supplementary Fig. 14c,d)*>*°. Moreover, depletion 
of PTEN in WT NSCs resulted in an elevated occupation 
of p-CREB and/or CBP as well as H3K4me3 at the PAX7 
promoter (Fig. 4k and Supplementary Fig. 14b). Using a 
GAL4-based reporter assay, we found that the co-activator 


activity of CBP was elevated in PTEN-deficient NSCs, which 
could be rescued by reintroduction of PTEN (Supplementary 
Fig. 14e). These results support a model that PTEN deficiency led 
to a de-repression of CREB/CBP machinery at the PAX7 
promoter, resulting in activation of PAX7_ transcription 
(Supplementary Fig. 14f). 

Next, we investigated whether over-activation of CREB/CBP 
upon PTEN deficiency could underlie the neoplastic trans- 
formation of NSC. For this purpose, we first repressed CBP/p300 
activity in PTEN-/- NSCs by shRNA. Silencing CBP/p300 
effectively decreased the migration of PTEN-deficient NSCs 
(Supplementary Fig. 14g,h). In addition, we compared the 
abilities of the WT, protein phosphatase-dead mutant (Y138L) 
and lipid phosphatase-dead PTEN mutant (G129E) of PTEN in 
rescuing the malignant phenotypes of PTEN-/- NSCs. 
Although complementation with either the WT or the G129E 
mutant PTEN effectively rescued the neoplastic phenotypes 
in PTEN-/- NSCs, mutation of PTEN at Y138L resulted in 
loss of rescue ability (Fig. 41), which is consistent with its inability 
to dephosphorylate CREB (Supplementary Fig. 14i). These 
observations indicated that the neoplasm-suppressing activity of 
PTEN in NSCs is at least in part dependent on its ability of trans- 
repressing CREB/CBP. 

Together, these data support a model that PTEN functions as a 
trans-repressor at the PAX7 gene promoter in NSCs; 
PTEN deficiency relieves the repression and activates CREB/ 
CBP-mediated transcription of PAX7 during the neoplastic 
transformation of NSCs. 


The PTEN-PAX7 axis regulates malignant features of GBM. 
Our results support the view that human NSCs are the potential 
cell-of-origin for GBM, especially in the context of PTEN 
deficiency. It is perceivable that GSCs, the tumorigenic popula- 
tion within GBM, inherit genetic trait from mutated NSCs that 
potentially mediates the malignancy of GBM. We established 
three lines of GSCs from primary patient GBM specimens based 
on the established methods?" These GSCs exhibited 
features characteristic of NSCs and could be cultured under 
neurosphere or monolayer conditions (Fig. 5a,b, Supplementary 
Fig. 15a,b and Supplementary table 1)°°. Western blotting and 
immunofluorescence showed that these GSCs expressed normal 
levels of PTEN protein (referred to as parental GSCs; 
Supplementary Fig. 15c). To create isogenic GSCs that are 
deficient for endogenous PTEN (referred to as GSCs-PTENÍ?"), 
the parental GSCs were transduced with a lentiviral vector 
encoding PTEN shRNA. Silencing PTEN in the parental GSCs 
resulted in a marked increase of migration capability (Fig. 5c 
and Supplementary Fig. 15d). Correspondingly, GSCs-PTEN!°Y 
exhibited increased levels of phosphorylated AKT, phos- 
phorylated CREB and PAX7 (Supplementary Fig. 15e,f). 
Knockdown of CBP/p300 in GSCs-PTENÍ?V not only impaired 
PAX7 expression but also compromised cellular migration 
potential (Supplementary Fig. 15g,h), whereas overexpression 
of PAX7 in the parental GSCs resulted in a more aggressive 
migratory phenotype (the resulting cells referred to as 
GSCs-PAX7hish), similar to that observed in GSCs-PTENIOW 
(Fig. 5d and Supplementary Fig. 15d,i). 

To establish a connection between PTEN and PAX7 in human 
GBMs, we examined the Cancer Genome Atlas (TCGA) data set 
of 532 GBMs^?, and found that GBMs with decreased PTEN 
expression demonstrated increased PAX7 transcript (Spearman 
correlation coefficient = — 0.28, P value = 1.0 x 107 !; Fig. 5e), 
therefore providing a strong evidence that PTEN expression 
negatively correlates with PAX7 transcripts in clinical 
GBM settings. 
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Figure 5 | PTEN-PAX7 axis functions in GSCs and GBMs. (a,b) Immunofluorescence analysis of the indicated antigens in GSCs (a) and their 
spontaneously differentiated derivatives (b). Scale bars, 1mm (phase) and 25 um (IF). (c,d) Cell migration assay indicating increased cellular mobility in 
PTEN-knocked down (c) or PAX7-overexpressed (d) GSC lines. n=3. ***P< 0.001 (t-test). (e) GBM microarray data were obtained from TCGA website 
and spearman correlation was calculated by all the samples which were plotted by red and blue dots. The gene expression level was log2 lowess 
normalized. PAX7 expression negatively correlated with PTEN (Spearman correlation coefficient — — 0.28, P value =1.0 x 107"). (fg) FACS (f) and 
TdT-mediated dUTP nick end labelling (g) analyses of cell apoptosis in MMC (7.5 uM, 24h) or vehicle-treated WT and PTEN~/~ NSCs. n= 6. *P« 0.05 (t- 
test); ***P < 0.001 (t-test); NS, not significant. (h) Representative bioluminescent images of PTEN -/ - NSCs implanted animals in the presence or absence 
of MMC treatment are shown at day 70 after injection (left). Right panel indicates average bioluminescent signals. n — 5 for dimethylsulphoxide (DMSO) 


treatment, and n— 8 for MMC treatment. ***P < 0.001 (t-test). 


Mitomycin C (MMO) triggers apoptosis in PTEN-deficient NSCs. 
Given that the PTEN deficiency is associated with increased 
metastasis, recurrence and drug resistance in GBM patients, 
discovery of compounds selectively killing PTEN-deficient GSCs 
may represent an important step forward towards precision 
therapy on GBM patients with known PTEN mutation(s). We 
therefore sought to utilize the isogenic NSC platform to screen 
compounds specifically kill PTEN-deficient NSCs while not the 
WT NSCs. By testing compounds with reported tumour-killing 
activity (Supplementary Table 2), we identified MMC, a clinically 
used anti-tumour agent with DNA-crosslinking activity ^^, 
as a potential candidate compound. MMC selectively killed 
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PTEN~'~ NSCs while showed minimal effect on WT NSCs, as 
evidenced by propidium iodide staining and TdT-mediated dUTP 
nick end labelling assays (Fig. 4£g and Supplementary Fig. 16a). 
Temozolomide, a clinically used anti-GBM drug, however, had a 
comparable effect on WT and PTEN- /- NSCs (Supplementary 
Fig. 16b). MMC was also found cytotoxic to NSCs overexpressing 
PAX7 (Supplementary Fig. 16c). 

To evaluate the efficacy of MMC in vivo, we performed 
intracranial administration of MMC to mice harbouring 
neoplasm derived from PTEN-/- NSC engrafts. MMC- 
treatment effectively repressed neoplastic growth from the 
PTENT’ NSC engrafts (Fig. 5h). 
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Discussion 
GBM is one of the most deadly forms of human cancers and has a 
median survival time of only ~ 14.6 months. Up to date, most 
studies of GBM pathogenesis rely on GBM cell lines, primary 
GSCs or GBM cells from fresh specimens, or mouse GBM 
models. These conventional source materials have undoubtedly 
aided in our understanding of GBM pathogenesis and facilitated 
development of novel therapies, however, have their drawbacks. 
For instance, GBM or GSC lines bear significant variations from 
line to line, such as marker expressions, differences in cellular 
behaviours in vitro or in vivo, and sensitivities to drugs. These 
differences can be attributed to different culture conditions, 
pathological stages or, more likely, from their distinct genetic 
backgrounds. The complex genetic mutations, including the 
primary, secondary or evolved ones in tumour suppressors and 
oncogenes, as well as their interactions constitute a roadblock for 
studying the oncogenic identity and cell-of-origin for GSCs, as 
well as for aetiological understanding of GBM. Although it has 
been well established that GSCs are directly linked to NSCs, it has 
been challenging to study how genetic mutations may contribute 
to malignant transformation of NSCs as well as initiation and 
progression of GBM. The primary GBM subtype is acute and of 
high-grade and frequently harbours loss-of-function mutations in 
PTEN gene. Mouse models have been employed to study the role 
of PTEN deficiency in NSCs1?995!, Given that species-difference 
exists between mouse and human, a human cellular platform to 
study GBM origin or pathogenesis associated with PTEN 
deficiency or other oncogenic events is thus urgently needed. 
It was not until very recently that two elegant human intestine 
tumour models generated using CRISPR/Cas9-based gene editing 
were reported?>°. Nevertheless, none of above-mentioned 
models investigated tumorigenesis in a genuine tumour- 
initiating niche. Our study for the first time uses genetically 
engineered human NSCs to study a specific genetic mutation on 
neoplastic transformation both in vitro and in intracranial 
tumour-initiating niches, and provided the proof-of-concept 
evidence that NSC is a most susceptible human adult stem cell 
type for PTEN-deficiency-induced oncogenic transformation. 
Nevertheless, it should be emphasized that GBM are frequently 
associated with more than one genetic mutations. Indeed, 
compared with the GSCs derived from clinically obtained GBM 
samples that potentially harbour multiple genetic mutations”, 
PTEN-deficient human NSCs show milder phenotypes in 
intracranial implantation assay, and do not show aggressive 
growth in non-neural microenvironment. We have shown that 
p53 insufficiency and PTEN deficiency resulted in a cooperative 
effect on NSC transformation, which is consistent with a recent 
mouse study indicating that simultaneous inactivation of Pten 
and p53 in mouse neural stem/progenitor cells led to malignant 
transformation and tumour formation!’. From this point of view, 
the isogenic NSC platform also provides a useful model to study 
other individual neoplastic mutations, or their combinations in 
mediating malignant human NSC transformation. 
Accumulating evidence indicated that PTEN plays a role in the 
nucleus through protein-protein interactions?^^^, We showed 
here that PTEN is a nucleus-localized protein in human NSCs. 
It associates with and inactivates CREB/CBP at PAX7 promoter, 
and mediates transcriptional repression of PAX7. To our 
knowledge, this is the first report to show that PTEN 
can trans-repress gene expression via the association with 
transcriptional activators on gene promoters. Notably, a recent 
report suggests that PTEN binds to Rad51 promoter and is 
required for its expression in MEFs’. In a separate study, PTEN 
is identified to be associated with histone H1 and accounts for 
maintenance of a condensed chromatin configuration**. Our 
findings, together with these observations, support the notion that 


nuclear PTEN may act as a key transcriptional regulator in a 
context-dependent manner. Given that PAX7 is among the 
earliest genes involved in neural progenitor commitment in the 
embryo, the elevated expression of PAX7 in PTEN-deficient GSCs 
may account for the low-differentiated histology of GBM, a 
common feature for more aggressive cancers. Together with the 
finding that PTEN is negatively correlated with PAX7 in clinical 
human GBM samples, it is likely that the overexpression of PAX7 
in GSCs is one of the key factors mediating GBM tumorigenesis 
or malignancy. 

Owing to its prevalence in GBM and other cancers, PTEN 
mutations have been regarded as a major target for drug discovery. 
Cell-based systems are advantageous for automate chemical 
screening or drug evaluation. Most of previous studies rely on 
cancer cells or primary cultures and use non-isogenic cells as 
control which are limited by confounding effect of different 
genetic backgrounds. When it comes to screening for specific drug 
to kill GSCs, human fetal NSCs from a different genetic 
background are the common cells to be used as a control". Our 
isogenic PTEN-knockout NSCs recapitulated certain features of 
GSCs, thus holding the potential to be used as a platform for 
screening ‘precision drugs’ against GSCs based on ‘PTEN 
deficiency. We found that MMC is a candidate compound, and 
preferentially targeting PTEN-deficient NSCs while sparing their 
PTEN-proficient counterparts. We further showed that 
transcranial injection of MMC could diminish neoplasia from 
PTEN-deficient NSCs. Nevertheless, although MMC is used in the 
clinic for other cancers, it does not cross the blood-brain barrier 
and has not been used for GBM. Thus, developing advanced 
MMC derivatives being capable of passing through blood-brain 
barrier and with reduced toxicity may represent a potential 
therapeutic strategy for GBM harbouring PTEN mutations, 
despite that the underlying pharmacological mechanism 
warrants further investigations. To our best knowledge, this is 
the first study using isogenic NSCs bearing a genetically 
engineered carcinogenic mutation as a potential drug discovery 
strategy. We propose here a simplified approach to examine ‘one 
mutation’ at a time, using human ESC (or iPSC)-derived NSCs 
with a well-known genetic background to screen or evaluate drugs 
being capable of targeting certain genetic mutations. 


Methods 

Antibodies and reagents. The following commercial primary antibodies were 
used: anti-PTEN (9559; Western blotting (WB),1:1,000; immunofluorescence (IF), 
1:200), anti-Phospho-PTEN (pSer380, 9551, 1:1,000), anti-Phospho-AKT (pSer473, 
4060, WB,1:1,000;IF, 1:200), anti-Phospho-AKT (pThr308, 2965, 1:1,000), anti- 
Phospho-GSK-3p (pSer9, 5558, 1:1,000), anti-Phospho-mTOR (pSer2448, 5536, 
1:1,000), anti-phospho-CREB (pSer133, 9198, 1:1,000), anti-CREB (9197, 1:1,000) 
antibodies were from Cell Signaling Technology; anti-NANOG (ab21624, 1:250), 
anti-Ki67 (ab16667, 1:1,000), anti-Nucleolin (ab22758, 1:200), anti-CBP (ab2832, 
1:1,000) antibodies were from Abcam; anti-OCT4 (sc-5279, 1:100), anti-SOX2 (sc- 
17320, 1:100), anti-Actin (sc-130301, 1:2,000), anti-GFAP (sc-9065, 1:1,000), anti- 
TRA-1-81 (sc-21706, 1:20) antibodies were from Santa Cruz Biotechnology; anti- 
Nestin (MAB5326,1:500), anti-Human nuclei (MAB1281, 1:1,000) antibodies was 
from Millipore; anti-MAP2 (m4403, 1:500), anti-Tujl (T2200, 1:500), anti-Flag 
(M2, 1:2,000), anti-SMA (A5228, 1:100), anti-AFP (A8452, 1:200), anti-tubulin 
(T5168, 1:2,000) antibodies were from Sigma; anti-CD34 (555822, 1:50), anti-CD43 
(560198, 1:50), anti-CD45 (555482, 1:50), anti-CD73 (550741, 1:50), anti-CD90 
(555595, 1:100) antibodies were from BD Biosciences; anti-CD105 (17-1057, 
1:100) antibody was from eBioscience; anti-PAX6 (PRB-278P, 1:500) antibody was 
from Covance; anti-PAX7 antibody (WB, 1:500; IF, 1:50) was from Atsushi 
Kawakami; anti-CD133 (130-080-801, 1:10) was from Miltenyi Biotec; MMC 
was purchased from Biomol International. Temozolomide was purchased 

from Sigma. 


Cell culture. H9 ESCs were purchased from WiCell Research Institute and cultured 
as previously described?^. Briefly, human fibroblast-derived iPSCs, H9 ESCs and 
their PTEN knockout derivatives were maintained on MMC inactivated MEFs in 
human ESC culture medium; whereas the feeder-free cultures were maintained on 
Matrigel- (BD Biosciences) coated plates in mTeSR medium (STEMCELL 
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Technologies). The human ESC culture medium contains: 80% DMEM/F-12 
(DMEM; Invitrogen), 20% knockout serum replacement (Invitrogen), 1% 
GlutaMAX (Invitrogen), 1% non-essential amino acids (Invitrogen), 1% Penicillin- 
Streptomycin (PS, Invitrogen), 55 «M B-mercaptoethanol and 10 ng ml ^ ! of human 
basic fibroblast growth factor (bFGF; Joint Protein Central). MSCs were grown on 
Gelatin-coated plates in culture medium containing 90% -MEM (Invitrogen), 1096 
fetal bovine serum (AusGeneX), 1% non-essential amino acid, 1% PS and 1 ng ml ^ 
of bFGF. NSCs were cultured in Neural Stem cell Maintenance Medium containing 
50% Advanced DMEM/F12, 50% Neurobasal, 1 x N2, 1 x B27, 2mM GlutaMAX, 
10 ng ml ^ ! of human leukemia inhibitory factor (hLIF), 3 uM CHIR99021 and 

2 uM SB431542. GSCs were grown in medium containing 90% Neurocult NSC Basal 
Medium (STEMCELL Technologies), 10% NeuroClut NSC Proliferation 
Supplements (STEMCELL Technologies), 20 ng ml ^ ! of EGF and 20 ngml ^! of 
bFGF and 2 ug ml ^ ! of heparin??. HEK293T cells were maintained in high-glucose 
DMEM containing 1096 fetal bovine serum (FBS, Gemini). No mycoplasma 
contamination was observed during cell culture. 


Generation of PTEN knockout (PTEN - / ~) human ESCs. TALENs against the 
human PTEN gene (PTEN-TALEN. Left and PTEN-TALEN Right) were pur- 
chased from Addgene (Plasmid #36764 Plasmid #36765). The donor plasmid for 
HR was described previously??. The following PCR primers were used for cloning 
the homology arms into the donor vector (PTEN-HR-Neo-pCR2.1): 

Left-F: 5'-CGGGGTACCGACCGTCCCTGCATTTCCCTCTACACT-3'; 

Left-R: 5'-CGCGGATCCCGGAATGGGGAGAAGACGAATAATCCTCC 

GAAC-3’; 

Right-F: 5'-ACGCGTCGACCCTGGTTGCAAGTGTCAA GCCACCGAT-3’; 

Right-R: 5‘-ATAAGAATGCGGCCGCCCCCATCCCTAATCAAAACCAATG 

TGTTGTACCTG-3’; 

H9 ESCs were pre-treated with ROCK inhibitor Y-27632 (Sigma) 24h before 
plasmid electroporation. On the day of electroporation, H9 ESCs were 
individualized with TrypLE (Invitrogen) and filtered through a 40-um cell strainer 
to remove cell clumps. For TALEN-mediated engineering of PTEN~/~ ESC lines, 
106 cells were counted and resuspended in 1 ml of MEF-conditioned medium 
containing 10 uM ROCK inhibitor Y-27632, then mixed with plasmids (8 ug 
PTEN-TALEN Left, 8 ug PTEN-TALEN_Right, 2 jig Neo-PTEN-HR-pCR2.1 
donor vector) and electroporated. After electroporation, cells were plated on MMC 
inactivated G418-resistent MEF cells. Two days later, cells were treated with 
50 ig ml - ! of G418 (Gibco) for selection. After 14-21 days of selection, 
G418-resistant clones were manually picked and transferred to 96-well plates and 
expanded for genotyping. Gene-targeting was determined by PCR-screening of 
genomic DNA from the drug-resistant clones with the following primers: 

P1: ,-AGGCTGTTACAGTCAAATCTCTGCGAACGAT-3^ 

P2: 5'-CCCCAAAGGCCTACCCGCTTCCATTGCTCA-3'; 

P3: 5'-CTACCTGCCCATTCGACCACCAAGCGA AACATC-3'; 

P4: 5-TGTCCACTGCCACAATTCGCATTACCAAACTCA-3'. 

Homozygous deletion of exon 1 of the PTEN gene was verified by a two-step 
PCR with the following primers: 

P5: 5-CCCAGACATGACAGCCATCATCAA AGAGATCG-3^; 

P6: 5'-AAAACGGTGAGCAGCAATTCTGGCATCACA-3'; 

P7: 5'-TTTGGGCCCCACACAGTGCAAGCTAAAGCAACCCCGTCT-3'; 

P8: 5'-CCGCTCGAGAACATCACTGGTCCTGGGCAAAATCGAGA-3'. 

The PTEN locus was further examined by Southern blotting after extracting 
genomic DNA from PCR-positive clones. The neomycin-resistance cassette was 
removed according to the previously described method?2??. Briefly, in order to 
remove the neomycin-resistance cassette, PTEN knock-out hESCs (5 x 106 cells) 
were electroporated with pCAG-Flpo-2A-puro vector (10 ug), and then plated on 
MMC-inactivated DR4-MEF cells. Two days later, puromycin (1 ug ml - t; 
Invitrogen) was added to the medium to enrich Flpo recombinase-expressed cells. 
The emerging colonies were expanded, and the removal of the neomycin-resistance 
cassette was verified by PCR using the following primers: P1: 5'-AGGCTGTTA 
CAGTCAAATCTCTGCGAACGAT-3' and P4: 5’-TGTCCACTGCCACAAT 
TCGCATTACCAAACTCA-3'. 


Flow cytometry analysis. For TRA-1-81 analysis, each 100 pl of cell suspension (10° 
cells) was used. Cells were incubated with the primary antibody mouse anti- TRA-1-81 
(Santa Cruz) and then the secondary Alexa Fluor 568-labelled Goat anti-Mouse IgM 
(Invitrogen). For CD133 analysis, each 100 ul of cell suspension (10° cells) was used. 
Cells were incubated with the primary antibody CD133/1(AC133)-PE (Miltenyi 
Biotec). For cell cycle analysis, cells were processed with the Click-iT EdU Flow 
Cytometry Assay Kits (Invitrogen) according to the manufacturer’s instructions. 
Briefly, the cells were harvested after incubation with EdU (10 uM) in medium. 100 pl 
of cell suspension (10° cells) for each sample were used. For the detection of cell 
apoptosis, each sample (10° cells) was subjected to Annexin V-PE staining??. Cells 
were examined by fluorescence-activated cell sorting (FACS) using a flow cytometer 
(FACSort; Becton, Dickinson and Company), and the cell cycle populations were 
determined using ModFit software (Verity Software House, Inc.). 


Generation and characterization of NSCs. Neural induction was based on pre- 
vious reports?022. In brief, ESCs cultured on MEF feeder cells were challenged with 


NIM-1 medium (50% Advanced DMEM/F12 (Invitrogen), 50% Neurobasal 
(Invitrogen), 1 x N2 (Invitrogen), 1 x B27 (Invitrogen), 2mM GlutaMAX 
(Invitrogen) and 10 ng ml ^ 1 of hLIF (Millipore), 4 uM CHIR99021 (Cellagentech), 
3 uM SB431542 (Cellagentech), 2 M Dorsomorphin (Sigma) and 0.1 uM 
Compound E (EMD Chemicals Inc.). Two days later, medium was switched to 
NIM-2 (50% Advanced DMEM/F12, 50% Neurobasal, 1 x N2, 1 x B27, 2 mM 
GlutaMAX, 10 ngml-! of hLIF, 441M CHIR99021, 3 uM SB431542 and 0.1 uM. 
Compound E) for additional 5 days. The cultures were then split onto Matrigel- 
coated plates and maintained in neural stem cell maintenance medium (NSMM) 
containing 50% Advanced DMEM/FI2, 5096 Neurobasal, 10 ng ml ^ ! of hLIF, 1 x 
N2, 1x B27, 2mM GlutaMAX, 3 uM CHIR99021 and 2 uM SB431542. 


Neuronal differentiation assay. 1 x 10° NSCs cultured on Matrigel-coated one 

well of six-well plate were switched to neuronal differentiation medium containing 
DMEM/F12, 1 x N2, 1x B27, 200 uM Ascorbic acid (Sigma), 400 uM dbcAMP 
(Sigma), 10 ng ml~! of GDNF (Peprotech) and 10 ngml-! of BDNF (Peprotech). 
Laminin (Sigma) was used to facilitate differentiation. Cells were cultured for 21 

days, and then characterized by MAP2 and Tujl immunostaining??. 


Bisulfite sequencing of the OCT4 promoter. Bisulfite conversion of DNA was 
carried out using CpGenome Fast DNA Modification Kit (Millipore) following the 
manufacturer's recommendation. About 1 ug of DNA was used as starting material. 
A genomic fragment of the OCT4 promoter was amplified using LA Taq Hot Start 
Version (TAKARA) as previously described”®. In brief, PCR products were purified 
by gel extraction using QIAquick columns (Qiagen), and subsequently cloned into 
the pEasy-T1 vector (Transgen). Thirteen clones from each sample were sequenced 
with the universal primer M13. 


Derivation of MSCs from ESCs. MSCs were differentiated from ESC lines as 
previously described”. In brief, groups of 10-14 EBs were plated on Matrigel- 
coated six-well plates in - MEM (Invitrogen) medium with 10% FBS (Hyclone), 
1% PS, 10 ng ml - ! of bFGF and 5 ng ml ~ ! of TGFB (HumanZyme). Cells were left 
for about 10 days until confluent MSC-like populations occurred. The resulting 
cells were passaged once and then analysed and sorted by FACS using various 
antibodies related to the MSC signatures. The cells after sorting were set to passage 
0 (PO). Expressions of MSC-specific markers were verified by immunostaining with 
antibodies against CD73 (BD), CD90 (BD) and CD105 (eBioscience). 


Isolation and culture of human GSCs. Three glioblastoma samples were obtained 
in the Beijing Tiantan Hospital, following approval from the institutional review 
board, and written informed consent was obtained from all patients. GSCs were 
isolated according to previously described methods with minor modifications*°*, 
Briefly, tumour tissues were dissociated into single cells with Accumax Cell 
Dissociation Solution (Millipore) for 30 min at 37°C. 2 x 10° single cells were initially 
seeded into one well of a six-well plate to allow formation of spheres/aggregates. In 
addition to sphere culture, GSCs were also cultured as monolayer in Matrigel-coated 
culture dishes. Cells were stably maintained either as spheres or as monolayer culture. 


Neurosphere assay. 2 x 104 single NSCs were seeded into one-well of an ultra- 
low attachment 96-well plate (Corning) to allow neurosphere formation. Medium 
was changed every two days. The diameters of neurospheres were measured using 
Image J. 


Luciferase reporter assays. A 981-bp genomic DNA fragment containing the 
PTEN docking site upstream of the TSS of PAX7 was cloned into the pGL3- 
Promoter vector (Promega). The DNA fragment was amplified by PCR using the 
following primers: 5'-CGGGGTACCCAAGTCTGGCCCAGCCCACACACTCA 
G-3' (forward primer); 5’-CCGCTCGAGGCCCTCCCCCGGGAGAGATTGG 
GAGATC-3' (reverse primer). Transient transfection with the luciferase reporter 
plasmids was performed in 12-well cell culture plates according to the manu- 
facturer’s protocol. Cells were lysed 48h post transfection, and firefly and Renilla 
luciferase activity were measured with a Synergy H1 Hybrid Reader(BioTek), and 
Renilla was used as an internal control. 


Immunofluorescence microscopy. Cells were fixed with 496 formaldehyde in PBS 
at room temperature (RT) for 20-30 min. After fixation, cells were permeabilized 
with 0.496 Triton X-100 in PBS for 5 min at RT. After incubation with 1096 FBS or 
BSA in PBS for 30 min, cells were incubated at RT for 1 h or at 4?C overnight with 
the primary antibody. Subsequently, incubation at RT for 20 min with the corre- 
sponding secondary antibody was carried out. Nuclei were stained with Hoechst 
33342 (Invitrogen). A Leica SP5 confocal was used for Immunofluorescence 
microscopy. 


Teratoma analysis. Teratoma assay was performed as described”. In brief, 
3 x 10$ ESCs were administrated subcutaneously into NOD/SCID mice (male, 6-8 
weeks). Mice were killed 6-12 weeks after injection, and teratomas were analysed 
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by immunofluorescence staining. All animal experiments were conducted with 
approval of the Institute of Biophysics, Chinese Academy of Science. 


SA-f-Gal staining. Staining was conducted as described previously”. In brief, 
cells were fixed at RT in 2% formaldehyde and 0.2% glutaraldehyde for 3 min. 
Fixed cells were then stained with fresh staining solution at 37°C overnight. The 
percentage of cells positive for SA-B-Gal signal were quantified and statistically 
analysed. 


Lentivirus production. Expression and purification of lentiviruses were performed 
as previously described with minor modifications”. For generation of lentiviral 
vectors encoding shRNA targeting PTEN, PAX7 or CBP/p300, corresponding 
shRNA oligos (Supplementary Table 3) were cloned into the pLVTHM vector 
pre-cleaved by Clal and Mlul. For lentiviral vectors expressing PTEN, PAX7 or 
luciferase (control), corresponding cDNAs were cloned into the pLE4 vector 

(a kind gift from Dr Tomoaki Hishida). For lentiviral packaging, HEK293T cells 
were co-transfected with lentiviral expression or shRNA vectors, as well as psPAX2 
(Addgene, 12260) and pMD2.G (Addgene, 12259). Lentiviral particles were 
collected after 48h and concentrated by ultracentrifugation at 19,400 g for 2.5h. 
Concentrated viruses were used for titre determination and cell transduction. 


Western blotting, IP and RT-qPCR. For Western blotting, cells were lysed in RIPA 
buffer (25 mM Tris-HCl pH 7.6, 150 mM NaCl, 1% NP-40, 1% sodium deoxycholate, 
0.1% SDS) with protease inhibitor cocktail (Roche). A BCA kit (Thermo Fisher 
Scientific) was used for protein quantification. About 20 ug of protein lysate was 
subjected to SDS-PAGE electrophoresis and electrotransferred to a PVDF 
membrane (Millipore). Then primary and HRP-conjugated secondary antibodies 
were used to visualize the target proteins. The quantification of western blot was 
performed with Image Lab software for ChemiDoc XRS system (Bio-Rad). For IP, 
~1x 107 cells was lysed on ice for 1h in NP40 lysis buffer (50 mM Tris-HCl, 
pH 8.0, 1% NP-40, 150 mM NaCl, 2mM EDTA, 1 mM PMSF, Complete EDTA-free 
Protease Inhibitor Cocktail). After adjusting protein concentrations to about 
4mgml ~ |, inputs were aliquoted, and lysates were incubated with anti-Flag agarose 
(Sigma) overnight at 4°C. After three washes in NP-40 lysis buffer, bound proteins 
were eluted with 1 x SDS loading buffer, and samples were analysed by western 
blotting with the indicated antibodies. For RT-qPCR, cellular total RNA was 
extracted by TRIzol (Invitrogen), and genomic DNA was removed by DNA-free Kit 
(Ambion), followed by cDNA synthesis with GoScript Reverse Transcription System 
(Promega). RT-qPCR was performed using iTaq Universal SYBR Green Supermix 
(Bio-Rad). Primer sequences are given in Supplementary Table 3. The uncropped 
images of the western blots and gels are shown in Supplementary Fig. 17. 


Cell viability assay. Cell viability was determined with CellTiter 96 AQueous 
One Solution Cell Proliferation Assay MTS (3-(4,5-dimethylthiazol-2-yl)-2,5- 
diphenyltetrazolium). 


In vitro cell migration assay. The cell migration of NSCs, MSCs and GSCs through 
the Matrigel (BD Bioscience) in vitro was performed as previously described??. First, 
the top chamber of the Transwells (24-well insert; pore size, 8 mm; Corning Costar) 
pre-coated with Matrigel was seeded with 2.5 x 10* cells in 100 ul basal medium 
(50% Advanced DMEM/F12, 50% Neurobasal for hNSCs, «MEM for hMSCs and 
GSC culture medium without EGF, FGF2 and Heparin for GSCs). The 
corresponding complete culture medium for each cell type was put into the bottom 
chambers. Cells that migrated through the Transwells were stained and counted. 


Single-cell clonal expansion assay. To determine the clonal expansion efficiency 
of NSCs and MSCs, 2,000 cells were seeded in a Matrigel-coated 12-well plate, and 
cultured until clear cell colonies formed. The relative colony number was then 
determined by crystal violet staining. 


Metabolomics evaluation. Metabolite profiling was carried out at Metabolon with 
LC-MS and GC-MS methods as described previously”. Briefly, 10° ESCs 

(PTEN +/+ versus PTEN-!-) and NSCs (PTEN */* versus PTEN-!-) were 
pelleted, and five replicates were stored at — 80?C. Samples were prepared and 
analysed by using Metabolon’s standard solvent extraction methods. 100 ul of each 
sample was extracted by methanol. After centrifugation, the supernatant consisting 
of metabolites was used for MS and LC-MS analyses. The extracts were equally split 
into two parts for analysis on LC/MS/MS (Thermo Fisher linear Ion Trap mass 
spectrometer with a Fourier Transform) and GC/MS (Mat-95 XP mass 
spectrometer) platforms. Data were analysed with Welch's two-sample t-test. 


Lactate measurements. Lactate secreted in the culture medium was collected and 
measured by a colorometric L-Lactate assay kit (Abcam) following the manu- 
facturer's instructions. 


Mouse brain in vivo imaging. Before intracranial implanting, NSCs were trans- 
duced with lentivirus expressing luciferase. 2 ul of a 1 x 107 cells per ml cell sus- 
pension in high-glucose PBS were injected to the hemi-striatum using stereotactic 
injection. The following coordinate parameters were used: antero-posterior — 0; 
medio-lateral = + 2.5 mm; dorso-ventral = — 3.5 mm. About 10 weeks post 
injection, mice were injected into peritoneal cavity with luciferin to trace tumour 
cells in vivo. The animals were next anaesthetized, and bioluminescence was 
imaged using the IVIS Lumina system (Perkin-Elmer). Images were taken every 
2 min for a course of 30 min (10s exposure per image). The mice brains were also 
imaged 1 and 2h post injection (10s exposure per image). Images were processed 
using Living Image software. Regions of interest were drawn based on each cell 
mass, and the total number of photons was recorded for each region of interest. 
Immunofluorescence staining and H&E staining were carried out on 5 um cryostat 
sections. 


Magnetic resonance imaging. MRI for mouse brain was performed on a Bruker's 
PharmaScan 7 T MRI imaging system. The system is equipped with standard 
actively shielded gradient system. A 72-mm birdcage resonator was used for 
excitation, and detection was done with a 30-mm surface coil. T2-weighted images 
were consecutively acquired by using a rapid-acquisition relaxation enhanced 
sequence. Animals were anaesthetized before MRI. 


RNA isolation and sequencing library construction. A total of 5-10 x 10? cells 
were used to isolate RNA with RNeasy Mini Kit (Qiagen). 1.5-2 ug of RNA was 
qualified by Fragment Analyzer (Advanced Analytical) and was used to construct 
sequencing libraries using the TruSeq RNA Sample Preparation Kit (Illumina) 
following the manufacture’s protocol. 


ChIP and sequencing. ChIP-seq was done following a previous protocol with 
minor modifications®. Briefly, cells were cross-linked by 1% (vol/vol) 
formaldehyde for 8 min at RT and then quenched by 125 mM glycine for 5 min at 
RT. Then, the cells were lysed in lysis buffer (50 mM Tris-HCl, 10 mM EDTA, 
1% SDS, pH8.0), which was followed by chromatin shearing using a Covaris S2 
instrument. The samples were then incubated overnight with Protein A-bead 
bound with the following antibodies separately: anti-PTEN (CST, 9559), anti-CBP 
(Abcam, ab2832), anti-phospho-CREB (CST, 9198) anti-H3K4me3 (CST, 9751s) 
and anti-H3K27me3 (CST, 9733). Normal rabbit IgG (Santa Cruz, SC-2027) or 
input was used as a negative control. To obtain the DNA of interest, chromatin 
with the target modification was eluted in elution buffer (20 mM Tris-HCl, 5mM 
EDTA, 50 mM NaCl, pH7.5), and then incubated with Proteinase K on a 
thermomixer at 1,300 r.p.m. for 2h at 68?C. Next, the DNA was extracted and 
purified via phenol-chloroform-isoamyl alcohol extraction and ethanol 
precipitation. Samples for sequencing including H3K4me3, H3K27me3 and IgG 
were constructed into libraries with NEBNext DNA Library Prep Reagent Set for 
Illumina (NEB) as described by the manufacturer. 


ChIP-PCR mapping of PTEN-associated genomic region. We mapped the 
PTEN-binding site at PAX7 promoter spanning from — 5,000 to + 3,251 with 
ChIP-PCR using different primer sets. In brief, to obtain the DNA region of interest, 
lysed samples were incubated overnight with Protein A-bead bound to anti-PTEN 
antibody. Then, extracted DNAs were subjected to PCR amplification of the different 
genomic regions by using different primer sets (Supplementary Table 3). 


Whole genome sequencing. Genomic DNA was extracted from 1 x 10° cells by 
DNeasy Blood & Tissue Kit (Qiagen) and sheared into 150-200 bp fragments by 
Covaris S2 and the quality was analysed via Fragment Analyzer (Advanced 
Analytical). The sequencing libraries were constructed using NEBNext DNA 
Library Prep Reagent Set for Illumina (NEB). 


RRBS and TA-RRBS. RRBS and TA-RRBS were performed as we previously 
reported®!, For RRBS, after digested with MspI, genomic DNA was size-selected 
and spiked in with 0.596 fully methylated lambda DNA. DNA was then treated with 
bisulfite after ligation with pre-methylated adapters and amplified with PfuTurbo 
Cx Hotstart DNA polymerase (Agilent Technologies). For TA-RRBS, after MspI 
digestion, 50-600 bp DNA fragments were first gel-selected, and then glucosylated 
with B-glucosyltransferase and oxidized by recombinant mouse Tet] (mTet1) 
subsequently. The final TA-RRBS libraries were generated after bisulfite treatment 
and amplification with PfuTurbo Cx Hotstart DNA polymerase as stated above. 


RNA-Seq data analysis. We aligned reads from RNA-Seq using TopHat® and 
Cufflinks? and differentially expressed genes were analysed as previously 
described®*, A total of 14,797, 14,782, 14,030 and 14,322 genes with FPKM >0.1 were 
detected in PTEN+/+ NSC, PTEN-/- NSC, PTEN^/* MSC and PTEN~!~ 
MSC, respectively. 281, 705 genes (q-value<0.05, P-value<0.05, FC>2) were 
significantly downregulated and 326, 1,139 genes (q-value<0.05, P-value<0.05, 
FC>2) were significantly upregulated in PTEN~'~ NSCs and MSCs, respectively. 
GO for enrichment of genes was assessed by DAVID®. 
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ChIP-seq data analysis. Sequencing reads from the ChIP-Seq experiment were 
mapped to the human reference genome (hg19) by Burrows-Wheeler Aligner 
(BWA)®, retaining only unique non-duplicate reads without mismatches 

in the first 15 bp. Peaks were then called by MACS by default parameters, 


using IgG as contro. 


167. Histone modification level was defined as non-duplicate 


reads in peaks (PR) normalized by total non-duplicate reads (NR), D = PR/NR. 
Modification level variation (PTEN -/- versus PTENt/+ >2) was assessed by 
the following command: macs2 bdgdiff -t1 k4wt.bdg -t2 k4ko.bdg -cl 
k4wt_control_lambda.bdg -c2 k4ko_control_lambda.bdg -d1 6191314 -d2 
6191314 -outdir./ — C 0.3 -g 200 -l 50. 


CREB-binding sites analysis. Promoter for CREB-binding sites analysis was 
defined as 3 kb upstream of TSS and 300 bp downstream of TSS. CREB-binding sites 
with P-value less than 1E-4 among the promoter regions were found by FIMO® 
using CREB motif downloaded from MEME database (MEME version 4.4). 


Genome-wide CNV analysis. Reads from whole-genome DNA sequencing were 
mapped to the human reference genome (hg19 from UCSC) by BWA. Signal of 

500 kb window CNV analysis was calculated by the following formula, 

C= 5*d*109/((500,000-£)*b), in which d is the base number in 500 kb that mapped 
to the reference, t is the number of N in the 500 kb window and b is the total bases 
that mapped to the reference. 


RRBS and TA-RRBS data analysis. First, sequencing reads were trimmed using 
trim galore v0.3.3 by the following command, trim galore --quality 20 --phred33 
--stringency 3 --gzip --length 50 --rrbs --paired --trim1 --output dir $out $raw1 
$raw2. Then trimmed reads were mapped to the hg19 reference genome by Bis- 


mark? and Bowtie 


7? using default parameters. Levels of 5hmC were determined by 


C/C + T of TA-RRBS and 5mC were determined by (C/(C + T) — 5hmC%). In our 
study, only genome regions with coverage > — 5 were taken into account. 


Statistical analysis. Data were analysed with Student's t-tests and were presented 
as mean + standard deviation (s.d. or s.e.m.). P «0.05 was considered as significant. 
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